Study of the induced potential produced by ultrashort pulses on metal surfaces 
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The influence of the induced potential on photoelectron emission from metal surfaces is studied 
for grazing incidence of ultra short laser pulses. To describe this process we introduce a distorted 
wave-method, the Surface Jellium-Volkov approach, which includes the perturbation on the emitted 
electron produced by both the laser and the induced fields. The method is applied to an Al(lll) 
surface contrasting the results with the numerical solution to the time-dependent Schrodinger equa- 
tion (TDSE). We found that SJV approach reproduces well the main features of emission spectra, 
accounting properly for effects originated by the induced potential. 

PACS numbers: 79.60.-i, 78.70.-g 



I. INTRODUCTION 

In the past few years developments in laser technology 
have made it possible to produce laser pulses with dura- 
tions in the sub-femtosecond scale [TJ |5J [31 [H |S] . This ad- 
vance in the experimental area opens up new branches in 
the research of the matter-radiation system [6| Ffl ISj [Qj ITO] . 
In particular, the investigation of photoelectron emission 
from surfaces due to incidence of short laser pulses gives 
the chance to understand a piece of the complicated puz- 
zle corresponding to electron dynamics at metal surfaces. 

In this article we investigate the photoelectron emis- 
sion produced when an ultrashort laser pulse impinges 
grazingly on a metal surface, focusing the attention on 
the role played by the surface induced potential. The 
induced potential is caused by the rearrangement of 
valence-band electrons due to the presence of the exter- 
nal electromagnetic field. This potential is expected not 
to affect appreciably electron emission for high frequen- 
cies of the laser pulse, for which surface electrons are 
not able to follow the fast fluctuations of the field. But 
for frequencies of the pulse close or lower than the sur- 
face plasmon frequency, the induced potential becomes 
comparable to the laser perturbation and its effect can- 
not be neglected. With this goal we introduce a simple 
model, named Surface Jellium- Volkov (SJV) approxima- 
tion, which includes information about the action of the 
surface induced potential, taking into account the main 
features of the process. 

The SJV approach is a time-dependent distorted wave 
method that makes use of the well-known Volkov phase 
[TT] to describe the interaction of the active electron with 
the laser and the induced fields, while the surface poten- 
tial is represented within the jellium model. This kind 
of one-active electron theories has been recently applied 
to study different laser-induced electron emission pro- 
cesses from metal surfaces, providing reasonable predic- 
tions [HI [131 [Ml [15]. To corroborate the vahdity of the 
proposed approximation, we compare SJV results with 



the numerical solution of the time-dependent Schrodinger 
equation (TDSE), in which the contribution of the sur- 
face induced potential is also included. The induced po- 
tential is here obtained from a linear response theory by 
considering a jellium model for a one-dimensional slab. 

With both methods - SJV and TDSE - we calculate the 
probability of electron emission from the valence band of 
an Al surface, considering different frequencies and du- 
rations of the laser pulse. We analyze in detail the ef- 
fect of the induced potential on electron distributions by 
comparing to values derived from the previous Impulsive 
Jellium- Volkov (UV) approximation |13j , which does not 
contain the dynamic response of the surface. 

The article is organized as follows. In Sec. II we 
present the theory, in Sec. Ill results are shown and 
discussed, and finally in Sec. IV conclusions are summa- 
rized. Atomic units are used throughout unless otherwise 
stated. 



II. THEORY 

Let us consider a laser pulse impinging grazingly on a 
metal surface (5*). As a consequence of the interaction, 
an electron (e) of the valence band of the solid, initially 
in the state i, is ejected above the vacuum level, ending 
in a final state /. The frame of reference is placed at 
the position of the crystal border, with the z axis per- 
pendicular to the surface, pointing towards the vacuum 
region. 

For this collision system we can write the Hamiltonian 
corresponding to the interacting electron as: 

H=Ho + Vl + V,nd, (1) 

where = — is the unperturbed Hamiltonian, 

with Vs the electron-surface potential, and Vl = r-F(t) 
represents the electron interaction with the laser field 
F(t) at the time t, expressed in the length gauge. In 
Eq.([l]), Vind denotes the surface induced potential, which 
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is originated by electronic density fluctuations produced 
by the external field. 

The electron interaction with the surface, V5, is here 
described with the jellium model, being Vs — —VcQ{—z) 
with Vc — Ep -\- Ew, where Ep is the Fermi energy, Ew 
is the work function and denotes the unitary Heaviside 
function. This simple surface model has proved to give an 
adequate description of the electron-surface interaction 
for electron excitations from the valence band of metal 
surfaces [13l [HI [15l HSl [17]. Within the jellium model 
the unperturbed electronic states, eigenstates of Hq, are 
written as: 
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(2) 



where the position vector of the active electron e is ex- 
pressed as r = (r^, z), with and z the components par- 
allel and perpendicular to the surface, respectively. The 
vector k = (k^,, kz) is the momentum measured inside the 
solid and i?k = k1/2 + e^^ corresponds to the electron en- 
ergy. The signs ± define the outgoing (-I-) and incoming 
(-) asymptotic conditions of the collision problem and the 
eigenfunctions 0^ (z) with eigenenergy e^^ are given in 
the Appendix of Ref. [18] . 

Taking into account the grazing incidence condition, 
together with the translational invariance of Vs in the 
direction parallel to the surface, we choose the electric 
field F {t) perpendicular to the surface plane, that is, 
along the z-axis. The temporal profile of the pulse is 
defined as: 



In this work we evaluate Tif by using two different 
methods: the SJV approximation and the numerical so- 
lution of the TDSE. Both of them are summarized below. 



A. Surface Jellium- Volkov approximation 

In the SJV theory, the final distorted state is rep- 
resented by the Surface Jellium Volkov wave function, 
which includes the actions of the laser field and the in- 
duced potential on the emitted electron, both described 
by means of the Volkov phase. The induced potential 
is derived from a linear response theory by using a one- 
dimensional jellium model [19]. It can be expressed as 
Vind{z,t) = z g{t) inside the solid, with the function g{t) 
numerically determined, while outside the solid - in the 
vacuum region- Vind{z, t) = 0. Hence, the final SJV wave 
function can be written as 

xf''^^^{v,t) = $^^(r,t)exp [iDl{kf,_,z,t)\ S,^nd{z,t), 

(5) 

where '^^^ is the unperturbed final state given by 

Eq.Q, which includes the asymptotic condition corre- 
sponding to emission towards the vacuum zone (external 
ionization process [TB]). In Eq. ([H]), the function DJ^ rep- 
resents the Volkov phase associated with the laser field, 
which is expressed as: 



D2{kf,,z,t) = - A-{t) ~ l3-{t)-kf, a-(t). (6) 



F (t) = Fq sm{ujt + ip) siu^ (wt/T) 



(3) 



for < i < T and elsewhere, where Fq is the maximum 
field strength, ut is the carrier frequency, ip represents the 
carrier envelope phase, and r determines the duration of 
the pulse. 

The differential probability of electron emission from 
the surface is expressed in terms of the transition matrix 
as: 



dP 
dk 



f 



dki Q{vp - ki) \Ti 
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where Tif is T-matrix element corresponding to the in- 
elastic transition kj^ k^ and k^ — (kj^,/c^j,) is the fi- 

/ 

takes into ac- 



nal electron momentum outside the solid, with /c', 
1/2 



In Eq. (41, ^ 2 



count the spin states and restricts the initial states 
to those contained within the Fermi sphere, with vp = 
{2EpY/'^. The angular distribution of emitted electron 
can be derived in a straightforward way from Eq. Q as 
d'^P/dEfd^f = k'f dP/dk'^, where Ef and flf are the 
final energy and solid angle, respectively, of the ejected 



electron and fc'. 



k', 



The temporal functions involved in Eq. (|6| are related 
to the vector potential A~ (t) , the ponderomotive energy 
P~{t) and the quiver amplitude a^{t) of the pulse, being 
defined as: 



A-{t) = -c /" dt'F{t'), 
J +00 



f3-{t) = (2c2)-i / dt'[A-{t')]\ 

J +00 

ft 



(7) 



a-{t) 



dt'A-{t'), 



-t-00 



with c the speed of light . In a similar way we express the 
function ^ind, which considers the action of the induced 
potential on the active electron, as 



cxp[i(z/c)A„^(i)] 



for z < 
for z > ' 
(8) 



with A^^j^ — — c J dt' g{t') the momentum transferred 

+00 

by the induced field. Note that the effect of the image 
charge of the emitted electron was not taken into account 
in the final distorted wave function X/^''^' because its 
contribution has been found negligible [5D]. 
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Employing the final SJV wave function given by Eq. 
5 1 within a time-dependent distorted- wave formalism 
the transition amplitude reads: 



where 



dt 



xf'^^-(i)|C/(t)|<i>+ 



it) 



(9) 



(10) 



represents the primary or collision (C) term, with 
U{z, t) = Vl{z, t) + Vind{z, t) the perturbation introduced 
by the laser and the induced fields and <I>j^ the unper- 
turbed initial state, given by Eq. ([2|. The second term 
of Eq. T^P(^\ is here called post-collision (PC) tran- 
sition amplitude, corresponding to the emission process 
after the pulse turns off at the time r. It reads: 



+00 



dt 



xf'''^-{t)\v,Ut)\KS^) 



(11) 



B. TDSE solution 

Replacing the semi-infinite jellium potential by the one 
corresponding to a one-dimensional slab of size a, Vgiab — 
—Vc 0(a/2 — z) Q{a/2 + z), and taking into account the 
symmetry of the system in the direction parallel to the 
surface, we can write the unperturbed eigenstates as 



*k.n(r,i) 
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-ipn{z)e 



(12) 



where now the functions ipn{z) are the discretized one- 
dimensional eigenstates of the slab potential. 

The time evolution of the electronic eigenstates un- 
der the laser pulse perturbation is governed by the one- 
dimensional time-dependent Schrodinger equation 

d 

l-Q^fniz, t) = H{z, t)Lpn{z, t) , (13) 

where the unperturbed part of the Hamiltonian H(z,t) 

is now Ho = -(l/2)(d2/dz2) + Vslab- 

The discrete time step evolution is given by the evolu- 
tion operator 



ipniz, t + At) = exp{—iAtH)ipn{z, t), 



(14) 



which is computed by using the Crank-Nicholson scheme, 
approximating the exponential by the Cayley form p2] 



exp{-iAtH) 



1 



iAt 



H 



1 



iAt 



H 



(15) 



This scheme is unitary, unconditionally stable, and accu- 
rate up to order (HAt)^. 



To obtain the transition amplitude we evolve every 
eigenstate within the Fermi sphere of the unperturbed 
slab, projecting then the evolved states over the dis- 
cretization box "continuum" states, (pj{z), 



t: 



{TDSE) _ 



{ip){z)\if, {z,t ^ oo)) . 



(16) 



Independence of the results with different slab sizes guar- 
antees that the used slab size accurately represents the 
semi-infinite medium. For the simulation box we have 
taken completely reflective walls as boundary conditions. 



III. RESULTS 

We applied the SJV and TDSE methods to study elec- 
tron emission from the valence band of an Al(lll) surface 
produced as a consequence of grazing incidence of ultra- 
short and intense laser pulses. As Aluminum is a typical 
metal surface, it will be considered as a benchmark for 
the theory. The Al(lll) is described by the following pa- 
rameters: the Fermi energy Ep = 0.414 a.u., the work 
function E\y = 0.156 a.u., and the surface plasmon fre- 
quency ujs — 0.4 a.u.. 

For the TDSE calculations a slab with a width of 
311.54 a.u. (142 Aluminum atomic layers), surrounded 
by 244.23 a.u. of vacuum on each side, was used. The 
grid sizes were Az = 0.1 a.u. for the spacial grid and 
At = 0.005 a.u. for the time grid. The time evolution 
was considered finished when the induced potential had 
decayed two orders of magnitude from its value at the 
moment the laser pulse was switched off, at t = t. The 
same criteria was used to evaluate the upper limit of the 
time integral of Eq. (111. Note that to compare the SJV 



and TDSE results it is necessary to take into account 
that the former theory includes the proper asymptotic 
conditions, distinguishing the external from the internal 
ionization processes, while the latter does not. Then, 
as a first estimation we weighted TDSE values with the 
fraction of electrons emitted towards the vacuum derived 
from the SJV model [13 . 

In this work we considered symmetric pulses, with 
ip = —LOT 12 -t- 7r/2 . The field strength was fixed as 
Fo = 0.001 a.u.(/ ~ 4 10^° W/cm^), which belongs to 
the perturbative regime, far from the saturation region 
and the damage threshold |23j i24j . In accord with results 
of a previous theory [T3] , the maximum of the emission 
probability corresponds to the angle 6'e = 90°, which co- 
incides with the orientation of the laser field. Therefore, 
all results presented here refer to this emission angle. 

Since the dynamic response of the surface is charac- 
terized by the surface plasmon frequency Wg, in order 
to investigate the influence of the induced potential we 
varied the carrier frequency w of the laser fleld around 
the value of cjj. We start considering laser pulses with 
several oscillations inside the envelope function, which 
correspond to the so-called multiphoton regime. In this 
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regime, related to a Keldysh parameter 7 = u}\/Ew / Fq 
[25j greater than the unity, the laser frequency tends to 
the photon energy and the electron spectrum displays 
maxima associated with the absorption of photons. 




0.5 

electron energy (a.u.) 
Figure 1 

FIG. 1: Differential electron emission probability, as a func- 
tion of the electron energy, for the emission angle 6^ — 90° . 
The parameters of the laser field are: Fq = 0.001 a.u., (a) 
uj = 0.7 a.u., T — 54 a.u., (b) uj = 0.4 a.u., r = 95 a.u. and 
(c) UJ = 0.2 a.u., r = 190 a.u. Solid (blue) line, SJV results; 
dashed (red) line, TDSE values, and dotted line, results of 
the IJV model. 



In Fig{l] six-cycle laser pulses with three different fre- 
quencies were considered: uj = 0.7, 0.4 and 0.2 a.u.. In all 
the cases, to analyze the effect of the surface response on 
the electronic spectra SJV and TDSE values were com- 
pared to data derived within the previous IJV approach 
|13j . which neglects the contribution of Vind- lu Fig{l] 
(a) we show the emission probability corresponding to 
the frequency uj = 0.7 a.u., which is higher than the 
surface plasmon frequency. For this frequency a good 
agreement between SJV and TDSE results is found. The 
SJV curve runs very close to TDSE values, showing only 
a small underestimation of TDSE results in high- velocity 
range. Note that both theories present a broad maxi- 
mum, which can be associated with the above threshold 
ionization process. From the comparison to values ob- 



tained within the IJV approximation, we observe that 
for this high frequency the induced potential produces 
only a slight increment of the probability at low electron 
energies, having small influence on the overall electronic 
spectrum. However, when uj becomes resonant with the 
surface plasmon frequency, as in Fig[l](b), the induced 
potential contributes greatly to increase the ionization 
probability in the whole energy range. Energy distribu- 
tions obtained with SJV and TDSE methods are more 
than one order of magnitude higher than the one derived 
from the IJV approach. In this case SJV results follows 
quantitatively well the behavior of the TDSE curve, de- 
scribing properly the positions of the multiphoton max- 
ima but underestimating TDSE probabilities around the 
second peak. Note that in this case Vind does not repre- 
sent a weak perturbation of the laser field, as shown in 
Figj2]^b) and it might be the origin of the observed dis- 
crepancy. In Fig|l] (c) we plot the emission probability 
for a laser field with a frequency uj — 0.2 a.u., lower than 
the plasmon one. Again, as in Fig[l] (a) SJV and TDSE 
results run very close to each other, displaying almost no 
differences with the IJV theory, which does not contain 
the induced potential. This indicates that the induced 
potential strongly affects emission spectra for frequen- 
cies resonant with the surface plasmon frequency, while 
for small deviations from this frequency it plays a minor 
role in the multiphoton ionization process. 

With the aim of examinating in detail the contribution 
of the induced potential, in Fig|2]we plot Vind as a func- 
tion of time, for a given position inside the solid and for 
the frequencies of FigJTJa) and (b). We observe that for 
u! — 0.7 a.u. the induced potential tends to follow the 
oscillations of the external field and its intensity steeply 
diminished when the pulse is turned off. Then, in this 
case the collective response of the medium produces only 
a weak effect on the electronic spectrum, as shown in 
FigjT|a). Whereas for laser frequencies near to ujs (Figj2] 
(b)) the process is dominated by the induced potential, 
which produces an increment of the emission probability, 
as observed in FigQb). 

Finally, in Fig|3] we study a six-cycle laser pulse with 
the frequency corresponding to the experimental value 
for the Ti:sapphire laser system [1] ( w = 0.057 a.u.). 
For this low frequency, almost one order of magnitude 
lower than the plasmon one, the surface response ap- 
proximates to the static limit and electronic fluctuations 
screen strongly the external field inside the solid. By 
comparing SJV and IJV results it is observed that in 
this case the induced potential contributes to reduce 
markedly the emission probability, up to two orders of 
magnitude at low electron energies. On the other hand, 
it should be noted that although the SJV theory describes 
properly the positions of multiphotonic maxima, it over- 
estimates the emission probability given by the TDSE 
method. Such a discrepancy, which arises when uj is 
lower than the mean energy of initial bound electrons, 
was also observed for other Volkov-type methods applied 
to photoionization of atomic targets pB] . 
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Figure 2 

FIG. 2: 2D-Representation of the induced potential, as a 
function of time and space. Laser pulse parameters are sim- 
ilar to Figjlj^a) and (b). Inset figures correspond to a given 
position inside the solid, with solid line, the laser pulse curve, 
and dashed line, the induced potential. 
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Figure 4 

FIG. 4: Similar to Fig[l] Half-cycle pulse with Fq = 0.001 
a.u., (a) frequency ui = 0.7 a.u. and duration r = 4.5 a.u., (b) 
frequency u — 0.2 a.u. and duration r = 16 a.u.. 




To complete the previous analysis we reduce the du- 
ration of the pulse in order to investigate the contribu- 
tion of the induced potential for photoelectron emission 



in the coUisional regime ||26j. In this regime, associated 
with half-cycle pulses, the electromagnetic field does not 
oscillate, producing a perturbation similar to the one 
resulting of the interaction with a swift ion impinging 
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grazingly on the surface (collision process). Notice that 
for such ultrashort pulses the carrier frequency to loses 
its meaning and the pulse can be characterized by the 
sudden momentum transferred to the ejected electron, 
Ap = -A-{0)/c ~ Fot/2 [37]. In Fig g] we plot elec- 
tron distributions for half-cycle pulses with two different 
durations r = 4.5 and 16 a.u. In both cases we found 
a good agreement between SJV and TDSE methods in 
the whole electron velocity range. Both theories present 
a pronounced maximum at low electron velocities, which 
does not appear in the electron distribution derived from 
the IJV approach, being produced by the induced po- 
tential. To understand the origin of this increment of 
the probability at low electron energies, in Figjsjwe plot 
again the induced potential for a given position inside the 
solid, now for the case of Figjl] (a). We observe that for 
half-cycle pulses, without oscillations, after the pulse has 
finished the induced potential still affects solid electrons 
during at least a hundred atomic units more. This effect 
is the main source of electrons emission at low velocities. 



A. Conclusions 

In the present work we have introduced the SJV ap- 
proximation, which allowed us to investigate the effects 
of the induced potential on the electron emission process. 
The proposed theory was compared to values derived 
from the numerical solution of the corresponding TDSE, 
displaying a good description of the main characteristics 



of photoemission spectra in the whole range of studied 
frequencies and durations of the laser pulse. From the 
comparison between SJV probabilities and those derived 
from the previous IJV approach, which does not include 
Vindi we conclude that the induced potential can play an 
important role in laser-induced electron emission from 
metal surfaces, as expected. For laser pulses with sev- 
eral oscillations inside the envelope, we found that the 
induced potential produces a considerable increment of 
the probability when the laser frequency is resonant with 
the surface plasmon one, but as to diminishes tending 
to the static case, the surface electronic density shields 
the laser field inside the solid, leading to a markedly 
reduction of the photoemission process. In addition, for 
electromagnetic pulses in the coUisional regime, the con- 
tribution of the surface induced potential after the pulse 
turns off gives rise to a maximum in the emission spec- 
trum at low energies. 
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